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Chapter 1 


Basic Concepts and Results in 
Aerodynamics 


1.1 Force Balance in Flight 


e The wings (and the horizontal tails in some cases) support the weight of the aircraft at subsonic 
speeds. 


e The rest of the aircraft just hangs from these “lifting surfaces.” 
e Of course the wings and tails themselves have weight. 


e On most aircraft, the wings contain most of the fuel. 


We can use Newton’s Laws of Motion to calculate the acceleration of an aircraft, and thus to decide 
how the forces on the aircraft must be balanced to make it go in a desired direction. 


Lift Distribution 


Wing Wt. 
Distribution 


Mg 


Figure 1.1: Force Balance in Flight 
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1.2 Straight And Level Steady Flight 


In Straight and Level Steady Flight, where all the accelerations are zero, Lift = Weight, and 
Thrust = Drag. 


L=Ww (1.1) 
p= 7 1.2) 
LIFT 
THRUST DRAG 
WEIGHT 


Figure 1.2: Straight and Level Steady Flight 


1.3 Basic Concepts and Results in Aerodynamics 


Freestream Vector: Velocity of the fluid far ahead of the object in the flow, undisturbed by the 
presence of the object. 


aircraft 


aircraft = 


Figure 1.3: Freestream Vector 
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For straight and level flight, 
L=W andT =D. 


Figure 1.4: Lift, weight, thrust and drag 


1.4 Planform Area S$ 


S is the area enclosed by a full-scale drawing of the outline of the wings when viewed from directly 
above or below. 


Figure 1.5: Planform Area S$ 


1.5 Dynamic Pressure in low-speed flow 


1 
dao = 5002, (1.3) 
p U |U_ |q q 
kg/m? | m/s | fps | N/m? | psf 
1.2 10 32.8 | 60 1.253 
1.2 
0.36 
Lift = (dynamic pressure) x (planform area) x (lift coefficient) 
L = QoSCL 
1 
L = 5 Uso SCL (1.6) 


The lift coefficient depends on how the lift was generated. Cz is usually around 0.1 to 1.4. 
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1.6 Lift Generation 


Generation of aerodynamic lift can be explained from Newton’s laws of motion. When flowing air 
is induced to turn, there is a rate of change of momentum along the direction perpendicular to 
the original direction. The net force which is required to cause this change, causes an equal and 
opposite reaction on the object which turns the flow. This reaction is the force acting on the object. 
Most simply there are three different “passive” methods to turn the flow. 


Rate of change of momentum = Force 


exerted on fluid 
=a <a> 
<== 2 
ee ed ————— 
<== 


Reaction exerted by fluid = Lift 


Figure 1.6: Lift Generation 


Irrotational Flow Around a Flat Plate -- ‘¥ Contours 
y 


rs As 7 


Figure 1.7: Airfoil At Angle of Attack: Streamline Pattern 
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15 degrees 


Display: = Streamlines 
http://www.windsofkansas.com/AZ2Fig5.jpg 


Figure 1.8: Flat Plate At Angle of Attack: Streamline Pattern 


How to turn the flow and generate lift: Camber 


——< 


Figure 1.9: How To Generate Lift Using Camber 
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1.7 Evolution of the Airfoil Shape 


Flat Plate, Stall - 


ae ~) Sy High drag too! 


Cambered Plate to Delay Stall 


oe 


Rounded leading edge gives larger range 
of angle of attack, but blunt trailing edge causes 
stall, and is ineffective in flow turning. 
Rounded leading edge to delay stall; 
sharp trailing edge to turn the attached flow. 


Se 


Figure 1.10: Evolution of the Airfoil Shape 


1.8 Vortex-Induced Lift 


———? 


Side Force 


Vortex 
Induced 
Lift 


Flow gets turned down by vortices. 
Roll Moment 


Figure 1.11: Vortex-Induced Lift 
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1.9 Vortex-Induced Lift and Delta Wings 


SP 


High-Aspect Ratio Wing elta Wing 


Vortex Lift 
Regime 


Figure 1.12: Vortex-Induced Lift and Delta Wings 


1.10 Airfoil (British: “aerofoil” ) 


“Airfoil” means “shape of a section of a wing.” It is a two-dimensional concept. 

Airfoils cannot fly: wings fly. 

Airfoil properties are used to calculate and design wing properties. 

Airfoil lift coefficient Cy, varies with angle of attack a. 

If the airfoil is cambered, the lift coefficient is positive even at zero angle of attack, and reaches zero 
only at some negative value of angle of attack: this is called the “zero-lift angle of attack,” ap. 

As the camber is increased, a9 becomes more negative. 

Thus airfoil lift coefficient is 


CL= Ta (7 29) (1.7) 


The lift-curve slope dC/da tells us how much lift increase we can expect to get from increasing 
angle of attack, from a given angle of attack. 
dC, /da < 27, where a is in radians. 
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1.11 Airfoils vs. Wings 


L' = Lift per unit span = qoocCy (1.8) 
D' = Drag per unit span = qoocCp (1.9) 


c: Airfoil chord, which is the line joining leading edge and trailing edge. 


Stall 


Note: Almost 
straight line 


O 


Figure 1.13: Cz curve 


1.12 Lift Curve Slope of an Airfoil 


Figure 1.14: Lift Curve Slope of an Airfoil 


“Thin-airfoil theory” in low-speed flows where Mach number is close to 0, proves that: 


d 
< < 27 in low-speed flows (1.10) 
ay 
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1.13. Exercise: Lift coefficient 


The angle of attack of an airfoil is 12 degrees. The lift curve slope is 5.8 per radian. Zero-lift angle 
of attack is -2 degrees. Find the lift coefficient. 


If the air density is 1/10 of sea-level standard, and the temperature is 20 °C higher than the 
standard sea-level, flight speed is 100 m/s and wing planform area is 30 m?, find the lift. 


Sea level standard temperature is 15 °C or 288 K. 


1.14 Pressure Coefficient 


The pressure coefficient is a way to express the pressure with respect to some reference pressure, as 
a “dimensionless” quantity. 


p=p p=p ins 
Cy = a ae | 1.11 
P™ 0.5pU2, is (=) ae 


Cy = 0 indicates the undisturbed freestream value of static pressure. 
Cp = 1 indicates a stagnation point. 
Cy < 0 indicates a suction region. 


Figure 1.15: Chordwise pressure distribution over an airfoil in low-speed flow 
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1.15 Flat Plate Streamlines and Pressure Distribution 


15 degrees 


Display: = Streamlines 


Output 
103.326 


Free Stream 


101.553 
q3.70 
49.007 


46.235 


Fig. 5 


http://www.windsofkansas.com/AZ2Fig5.jpg 


Figure 1.16: Flat Plate Streamlines and Pressure Distribution 


1.16 Exercise: Pressure Coefficient 
e What is the pressure coefficient at the stagnation point of an airfoil section? 
e What is the pressure coefficient on a flat surface aligned with the freestream? 


e Cy at the suction peak of an airfoil is -1.2. What is the pressure there as a percentage of the 
freestream static pressure? 


e What is the velocity at this point, as a percentage of the freestream velocity? 
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1.17 Drag 


Ideally, nothing happens to the fluid when an object goes through, except that it gets moved aside, 
or maybe turned a bit, but then recovers. 

In practice, some effect always remains. Some flow gets pulled along the direction of the object 
until its kinetic energy is dissipated as heat. Or there is a region of swirling flow left behind, whose 
energy is also dissipated as heat. The force resulting from these things is drag. It acts along the 
freestream direction. 

The drag is given by: 


D = qxnSCp (1.12) 
Lift to Drag Ratio: 
CL 
L/D=— 1.13 
[D == (1.13) 


We want our airplanes to have as high L/D as possible! 


> - x 


Figure 1.17: Swirling flow left behind the airfoil 


The profile drag of an airfoil of chord 1 unit is about the same as that of a circular cylinder 
whose diameter is only 0.005 unit. 


oe 


Airfoil of chord c 


U ; Cylinder of diameter 0.005c 


Figure 1.18: Profile Drag of a Cylinder vs. Airfoil 


Streamlining: At high Reynolds numbers, the drag of an airfoil with chord c is only as much as that 
of a cylinder of diameter 0.005c. From Covert, Eugene C., “The Legacy of the Wright Brothers and 
Its Future,” 1997 Wright Brothers Lecture. 
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But we want our space reentry capsules and parachutes to have very low L/D. 


http://ares.jsc.nasa.gov/HumanExplore/Exploration/EXLibrary /images/Pics/MarsScenariol/09Chute.gif 


http://serc.carleton.edu/images/NAGTWorkshops/careerprep/jobsearch/Parachute_300.jpg 


Figure 1.19: Examples of low L/D 


1.18 Drag Coefficient 


The drag is given by D = 0.5pU2,SCp. 
The drag coefficient in low-speed flow is composed of 3 parts: 


Co = Coo + Cpfriction + Cdi (1.14) 


when Cpo is the parasite drag, which is independent of lift. It is usually due to the losses of stag- 
nation pressure which occur when part of the flow separates somewhere along the wing or body 
surface. In high speed flight, the effect of shocks and wave drag must be added to this, and becomes 
dominant source of drag. 


Example: 

Cpo of a small airliner is 0.018. Wing aspect ratio is 6. Assume spanwise efficiency is 1.0. Lift 
coefficient is 0.5. Find the total drag coefficient. 

If the density is 1 kg/m? and speed is 200 m/s, find the drag. 


1.19 Drag Calculation 


Drag = (Dynamic Pressure) x (Reference Area) x (Drag Coefficient) (1.15) 
D = q0SCp (1.16) 
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D 
rs = SCp is called “Equivalent Flat Plate Area” (1.17) 


Note: Unlike lift, which people associate with the wings and hence planform area S', everything 
produces drag. So when the drag of each component is calculated, the drag coefficient of that is 
based on the “referenced area” for that component. For instance, for a fuselage the reference area 
is the cross-section area. 

Thus when everything is summed up, you use the concept of “Equivalent Flat Plate Area” rather 
than Cp to make sure we don’t use different reference areas. 


1.20 Example: Equivalent Flat Plate Area 


Subscript “0” indicates “zero-lift” value. 


Component Cpo | Area Do/q, ft? 
Fuselage + Engine Nacelles | 0.03 | C/S area: 64 ft? ? 
Wings 0.003 | Planform area: 120 ft? ? 
Tails 0.004 | Tail planform area: 15 ft? ? 
Landing Gear 12 Strut diameterx Length = 2/12x5 | ? 


Total Do/gq [ft?] of the airplane = ? 
C'po based on planform area = ? 


Let’s say that the airplane is flying at minimum drag: Lift-induced drag = Profile drag. Let’s 
say that the aircraft weight is 5000 Ibs. What is the L/D? 


1.21 Lift-Induced Drag and Aspect Ratio 


At the ends of the wings, the pressure difference between the upper and lower sides is lost, as the 
flow rolls up into a vortex. 


Inboard vortex sheet 


Tip-vortices 


Figure 1.20: Tip vortices 
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The aspect ratio of a wing is defined as: 


b2 
AR => (1.18) 


where b is the wing span and S is the wing planform area. 


Span b 


| 


Figure 1.21: Wing span b 


1.22 Effects of Finite Aspect Ratio 


1. The overall lift is reduced, relative to the airfoil lift value predicted for a section of an infinite 
wing. 


2. The lift vector is tilted back, so that an “induced drag” is created. 


Drag, D 


itt, L Normal Force, F 


ind 


Figure 1.22: Induced drag 


Both of these (usually undesirable) effects are reduced by increasing the aspect ratio of the wing. 


1.23 Points to Note on Finite Wings 


1. Since induced drag is directly related to the lift it can be calculated by the same mathematical 
formulation used to calculate lift. 


2. Does not require consideration of viscosity. 
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3. No induced drag on 2-D airfoils under steady conditions. 


4. At finite Aspect ratios, 
Ch 


Cbi = TAR)e 


(1.19) 


5. Ideal elliptic lift distribution implies minimum induced drag, i.e. spanwise efficiency e = 1. 


Note: In the 2-D limit (airfoil) there is no lift-induced drag in incompressible flow. There is always 
“profile” drag due to viscous effects (friction drag) and flow separation (pressure drag). 


1.24 Effect of Aspect Ratio on Wing Lift Curve Slope 


=a= (1.20) 
da a 
s Er vane 
dual aay 
ie “1 (1.22) 


1.25 Typical L/D values of different types of aircraft 
e Transport aircraft: (L/D) 16 - 20 
e Fighters (L/D) cruise 10 - 16 
© Supersonic transport: (L/D) 11 


e Hypersonic aircraft: (L/D) 1 (lift is a minor problem here) 


1.26 Example on L/D and Wing Loading 


Example: 

An aircraft has a wing loading (W/S) of 130 pounds per square foot (6233 N/m?), aspect ratio 
of 7.667, and wing span of 60.96 m. We'll assume that its spanwise efficiency factor will be 0.99. 
Let’s assume that the profile drag coefficient is given by Cpg = 0.015. Thus, for maximum Lift-to- 
Drag ratio (minimum drag, and the lift is always equal to the weight for straight and level flight), 
Cp; = Cpo = 0.015. 

The corresponding Cy, is calculated as 0.598, and the dynamic pressure is 10423 N/m?. At 
11,000 meters in the Standard Atmosphere, density is 0.36 kg/m, so that the flight speed is 240.64 
m/s. 

Note: In practice, the Cpo9 might change with flight Mach number, for high-speed flight. This 
is not taken into account in the above. 
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1.27 Aerodynamics Summary 


Lift is force perpendicular to the flow direction, due to 
pressure differences across surfaces. 3 ways of generating 
lift are: 


a) angle of attack 
b) camber 
c) vortex-induced lift 


An infinite (2-dimensional) wing is entirely described 
by its airfoil section. Finite wings have less lift than 
corresponding span-lengths of an infinite wing at the 
same angle of attack, and also have lift-induced drag. 
The total drag is composed of profile drag, which does 
not vary with lift, and induced drag, which rises as the 
square of the lift coefficient. 

To fly an airplane of a given weight, straight and level, 
the condition for minimum drag (maximum lift-to-drag 
ratio) is that the profile drag coefficient is the same as 
the induced drag coefficient. 


1.28 Speed for Minimum Drag 


Total drag is composed of a part which depends on lift, and one that does not. 


D = Do + Dj = (Cpo + Cni) (0.5pU2,S) (1.23) 
Cr 2 
D= —_ 1.24 
Let us consider what it takes to keep L = W. 
W =L=q~SCr (1.25) 
WwW 
= —_ 1.26 
me dooS oe 
So F 
W 1 S 
D = qo 1.2 
estoo+ (5) wana oe 
dD W\? 1 S$ 
dD _ = 1 
ina = 8° - (5) ames ~° 


i.e. Cpo = Cp; or, Minimum Total Drag = Twice Zero-lift Drag 


... This is a remarkable result. It means that: 
AIRCRAFT, UNLIKE OTHER FORMS OF TRANSPORTATION, HAVE A DEFINITE SPEED 
FOR MINIMUM DRAG! 

To fly an airplane of a given weight, straight and level, the condition for minimum drag (maxi- 
mum lift-to-drag ratio) is that the profile drag coefficient is the same as the induced drag coefficient. 
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Speed for minimum drag 


Flight Speed U 


Figure 1.23: Profile Drag, Induced Drag and Total Drag versus Speed 


1.29 Kutta-Jowkowsky Theorem 


Derived from Newton’s 2nd law. 
Derived by observing that there are 2 ways that lift can be explained: 


1. Due to circulation (flow turning) 
Lift per unit span L’ = pUx x I. Here [ is the Bound Circulation. 


I 


Figure 1.24: Kutta-Jowkowsky Theorem 


2. Due to compression / expansion of the flow (density changes due to velocity changes): 
This cannot happen in incompressible flow, which is what we study in this course. 


= 


Figure 1.25: Compression and expansion of the flow 
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Chapter 2 


Introduction to Fluid Mechanics 


2.1 Types of Fluid Motion 


Fluids, like most other forms of matter, are made up of tiny particles (molecules), which are sepa- 
rated by large spaces. 


> ss 
.» * 


Figure 2.1: Molecular Kinetics 


How many molecules are there in a cubic centimeter? 


Avogadro’s # is 6.022115 x 1073. 
Gives # of atoms (or molecules) per gram-mole of the gas. : 

1 gram-mole of air is roughly 28.97 grams. Molecular weight M is 28.97 grams/mole, or kg/kg- 
mole. 


18 
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Density of air (mass/volume) can be found from the Perfect Gas Law. Let’s say, we have air at 
sea level standard conditions. Pressure is 1.01325 x 10° N/m?, and temperature is 288 K. 


== (2.1) 
pS RE 
PM 
RT 
101325 x 28.97 
8313 x 288 


Density is 1.225 kg/m?. 
So how many grams in a cubic centimeter? 1.225 x 107°. 


1.225 x 10-3 
How many molecules? acne x 6.0221415 x 1078 = 2.55 x 107°. 


So we can treat fluid as a “continuum”: 
a medium which is pretty uniform in properties at the smallest scales of interest to us. 


Fluid Motion 


In the rest of this course we will treat fluids as “continuous media,” and forget about those molecules 
zipping about at random. 

The smallest “packet” of fluid which we will now consider will have many millions of molecules. 
The motion of these packets is the net motion of all those molecules zipping around within the 
packet (think of a school bus moving at 5 mph with thirty middle-schoolers conducting a football 
game inside using each other’s backpacks. ) 


Streamline 


A streamline is a curve whose tangent at any point is in the direction of the velocity vector at that 
point. 
NOTE: if the flow is not steady, the streamline has little use or significance. 


Figure 2.2: Streamline 
If ds’ is an elemental vector along the streamline, 
dz = dai + dyj + dzk (2.4) 


U =uitvj+wk (2.5) 
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7 1 oj ok 
dsxU=|dx dy dz (2.6) 
“ ov w 
=> i(wdy — vdz) + j(udz — wdx) + k(vdx — udy) = 0 (2:7) 
Equation to a streamlines is: 
dz ow 
= ae ee 2: 
wdy = vdz a (2.8) 
udz = wdz ee (2.9) 
7 dx u : 


Four Basic Types of Fluid Motion 


As a “packet” of fluid moves along, a combination of 4 things can happen to it: translation, 
dilatation, rotation and shear. Thus any fluid motion can be described as some combination of 
these. 


1) Translation: motion of the center of mass 
This is characterized by the velocity 


U =uitvuj+wk (2.10) 
P 

= — 2.11 

P= pr (2.11) 


Figure 2.3: Translation 


2) Dilatation: volume change 
ws ea (2.12) 


Also called “Divergence of the velocity vector.” 


Figure 2.4: Dilatation 
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3) Rotation: about one, two or three axes 


cl dw dv\-+ dus dw\ + dv du\ >= 
Fs ; ) | 2.1 
w= 9 (+ ak (+ yi (= =) 7 ar) 


vorticity ¢ = 2uw =V x U (2.14) 


Note that rotation is a vector. We usually use a quantity called “vorticity” ¢ which is twice the 
rotation vector, to describe the amount of rotation in flows. 


Figure 2.5: Rotation 


4) Shear Strain 
Strain is defined as a change in length per unit length. Rate of strain is thus given by change 
in velocity per unit length. Shear exists when there is a gradient of velocity along the direction 
perpendicular to streamline. 
Ov Ou Ow Ov Ou Ow 
€ = — + _ € = — + — € = — + — 2.15 
"Ox Oy ve Oy Oz “Oz Ox ( ) 
The quantity “strain rate” has nine components. They are: 
Exes Eyys €zz These are “normal strain rate” components. 
Exys €yz, €za These are “shear strain rate” components. 
Eyas Ezy, €xz These are equal to the corresponding “shear strain rate” components above. 


2.2 Circulation [= — $. U- dl 


Defined as an integral around a closed contour “c.” The negative sign is included such that positive 
circulation on a body corresponds to positive lift, and the integral is evaluated counter-clockwise. 
From the preceding discussion, we see that will be zero unless there is some vorticity, contained 
within the contour. 

T is an extremely useful quantity: it helps us calculate lift, vortex strength, etc. 


Important points: 
The circulation around a closed contour with net rotation and/or shear will be non-zero. It is, how- 
ever, always possible to have a combination of rotation and/or shear that gives a zero circulation. 


Note: In the above, the definition is based on a line integral around a contour, which is in a 
given plane. This works when we know the plane of the rotation, and that the flow is entirely 
2-dimensional. To deal with general problems, we can consider I to be a vector aligned along the 
axis of rotation. 
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Figure 2.6: Circulation 


2.3. Kutta-Jowkowsky Theorem (revisited) 


Derived from Newton’s 2nd law. 
Derived by observing that there are 2 ways that lift can be explained: 


1. Due to circulation (flow turning) 
Lift per unit span L’/ = pU, x I. Here [ is the Bound Circulation. 


I 


Figure 2.7: Kutta-Jowkowsky Theorem 


2. Due to compression / expansion of the flow (density changes due to velocity changes): 
This cannot happen in incompressible flow, which is what we study in this course. 


——= 


Figure 2.8: Compression and expansion of the flow 


Chapter 3 


Conservation Equations of Fluid 
Mechanics 


3.1 Conservation of Mass 


If we have a fluid going in and coming out of a given region of interest, (a “control volume”), we 
can say for sure that (what goes in per unit time) = (what comes out per unit time) + (what 
accumulates inside per unit time). 


Figure 3.1: Conservation of Mass 


; ; ; AM 
my =mMm2+m3+ | —— (al) 
At /inside 
(Mass going in per second) = (Sum of masses going out per second) 
+(Mass accumulated inside in 1 second) (3.2) 


In general, mass may be going in and/or out everywhere across an (imaginary?) surface enclosing 
the space in which you are interested. Also, the velocity of the inflow/outflow may be nonuniform, 
and in some arbitrary direction. 


Mass flowing out of the surface of the control volume per unit time = fp pla. ads (3.3) 
S 


where dS is a small element of the total surface area S. An “integral” is a neat way of saying, 
“collect all the little bits and add them up.” 
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Law of conservation of mass for a control volume becomes 


ff oa-mas + fff pav =0 (3.4) 
Ss Vv 


“Conservation” of anything can be expressed this way, as we will see. Here the r.h.s. is zero: mass 
can’t get converted to anything. 


3.2 Example: Closed circuit wind tunnel 


Figure 3.2: Closed-Circuit Low Speed Wind Tunnel 


Settling Chamber velocity = 15 fps, uniform; steady. Chamber diameter: 20 feet. Test section 
is a 9-foot diameter duct. Our “control volume” has one face in the settling chamber, other face 
in the test section. Assume that no flow can escape out the sides of the control volume. “Steady” 
means that at any given point in the flow, the properties don’t change as time changes. 


6) 

3y (anything) = (3.5) 
The “integral form of the mass conservation equation” becomes: 

fp p(a- ad =0 (3.6) 

S 


Subscript “1” refers to settling chamber face of control volume, subscript “2” refers to test section. 
Air cannot escape out of the sides of the “control volume,” so the integral reduces to just 


—p, AU, = p2A2U2 = 0 (3.7) 


Negative sign: Look at direction of vector normal to each face, with respect to the direction of U. 


Thus 
_ mA 


p2A2 


Since the density change can be neglected, and each area is a circle, 


Uy (3.8) 


50 \7 
Us = (>) x 15 = 74.07 ft/s (3.9) 
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3.3. Conservation of Momentum 


Newton’s 2nd law of Motion: “Rate of Change of Momentum = Force” 
or “Net force acting on a system = time rate of change of momentum of system” 
Note units: 


Momentum = Mass x Velocity (3.10) 

= Density x Volume x Velocity (3.11) 

= Density x Area x Distance x Velocity (3,12) 
Momentumperunittime = Density x Velocity x Velocity x Area (3.13) 


What is being conserved? “momentum per unit volume”. So the conservation equation has 
2 terms on the left hand side: 


1. “rate of change of the quantity, integrated over the whole control volume”. This is the “un- 
steady term”. 


2. “net outflow per unit time of the quantity across the surfaces of the control volume, integrated 
over the whole control surface”. This is the convection term. 


On the right hand side, we have surface and volume integrals of the quantities to which our 
“conserved quantity” may have changed. In this case, when momentum disappears, the rate of 
change of momentum produces forces of various kinds. 


[J (oay(a- aya A + a [[femav = fff force (3.14) 
C.S. C.V. 


vol 


3.4 Body and Surface Forces 


a) Body forces 


Gravitational force per unit volume = pg (3.15) 


(Important for liquid flows, and buoyant gas flows). 

Electromagnetic force: it would be something like p(B x €). 

Important in ion propulsion, spark plugs, plasmas. Generally, we neglect body forces in aerody- 
namics. 


ay 


Surface forces 

Pressure is the result of molecules, flying about at random, crossing (or colliding with) the 
surface, thus transferring momentum. So it acts normal to the surface. It must be related to # 
of molecules/volume (thus related to density), and speed of random motion of molecules (thus 
related to temperature). Thus, pressure is related to the product of density and temperature. 


This relation is a simple proportionality (multiply by a constant) where R is the gas constant. 
This is called the “perfect gas law.” 


So force due to pressure is — ffs (pri)dS because pressure acts on the surface, opposite to 7. 
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3.5 Shear Stress vs. Rate-of-Strain Relations 


tan’() 
stress 
stress 
d ; 
strain oo 
Figure 3.3: Stress and Strain 
Ou Ov 
ay as PA 3.16 


a Oz Oy 


(3.17) 


Ou Ow 


The direction is indicated by the second subscript. The units of stress are Force/Area, such as 
N/m’, psi, or psf. 

These can be written in compact form. Here u is used to represent any of the velocity com- 
ponents, and x is used to represent any of the spatial coordinates. i and 7 representing x,y,z in 


turn. 
ue Ou; ft Ou; 
9 B\ Oa; | Oa; 


The normal stress, on the other hand, is simply the pressure. 


(3.19) 


3.6 Normal Stress 


Here “normal” is a fashionable term for “perpendicular to a surface,” as opposed to “tangential,” 
not as opposed to “abnormal”... 
The normal stress is simply due to the pressure. 


Orr = —pi (3.20) 
Oy = —Pi (3.21) 
O22 = —pk (3.22) 


The unit vectors are included to show the direction of each of these stress components. 
Don’t fluids resist the strain rate due to compression and expansion? 
They do, and a relation between normal stress and rate of normal strain can be written, just like 
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the shear strain relations on the previous slide, and added on to the above expressions. However, 
in gases, the coefficient of resistance to normal rate of strain is very small, so that this “bulk 
viscosity” does not become significant unless the strain rate is extremely high. 

This occurs, for instance, across a shock wave, where velocity decreases by a large factor in the 
space needed for a few (under 10) collisions between molecules. 

In most other situations, it is a usual assumption that “Bulk Viscosity is Zero” in a gas so 
that the normal stress is attributed to just the pressure. 


3.7 Energy Equation 


In addition to mass and momentum, we can use the fact that “energy is neither created nor 
destroyed, but can change form.” 

The “conserved quantity” is energy per unit volume. The lhs has the terms describing time rate 
of change, and “flux” across the surfaces of the control volume. The rhs has the terms describing 
what changes the energy in the control volume: work and heat transfer. We use the first law of 
thermodynamics, which says that if you do work or release heat, your energy is exhausted. We will 
leave detailed study of this equation to the course on compressible flow and thermodynamics. The 
measure of total energy in a flow is the “stagnation enthalpy,” which is related to the stagnation 
value of temperature. So the energy equation reduces to: “Stagnation temperature goes up if 
work or heat are added to the flowing fluid.” 

In steady flow, with no work being added, and no heat being added (adiabatic flow), the rhs is 
zero, so that the energy equation reduces to “Stagnation Temperature is Constant”: 


uz 


To =T+H+ = constant (3.23) 


2Cp = 
The c, here is not pressure coefficient, but “specific heat” at constant pressure. 


ky 
C= 3.24 
pay (3.24) 
For diatomic gases such as air, at usual temperatures that we encounter in low speed aerodynamics, 
the “gamma” has the value 1.4. Thus c, = 3.5R. 


3.8 Converting the Integral Form to the Differential Form of the 
Conservation Equations 


In this section we will first convert the integral conservation equations to a form suitable to apply 
at each point, so that we can track changes from one point to another. 

So we have to go from the “integral form” over a control volume, to a “differential form” which 
deals with small changes from point to point. 
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of 


oe 
i. wN : a 
Joy 
—————_ cp . Ff 
dt ? OX 


<S 


Figure 3.4: Integral over a region vs. gradient around a point 


3.9 Relations between integrals over lines, surfaces and volumes 


Why? Because we would like to bring everything under the same integral sign, set the whole integral 
to zero, and relate the quantities inside, thus getting rid of all the integrals: 


dS dV dl 3.25 
ffows [fo £0 (3.25) 


Stokes’ Theorem 
Relates line integral over a closed contour to surface integral over the surface enclosed by the contour 


faad- pv x @) a5 (3.26) 


S 


This can also be written with the area element dS as a scalar, and the unit normal vector indicating 
its orientation. 


fa di = {pw x i) ads (3.27) 


S 


Divergence Theorem 


Relates surface integral over a control surface to the volume integral over the volume enclosed by 


the control surface. 
ff (z-a)dS = {fv -i)dV (3.28) 
S V 


Gradient Theorem 


Relates surface integral of a scalar a control surface to the volume integral of the gradient of the 
scalar over the volume enclosed by the control surface. 


fp pd§ = {[fcvpyav (3.29) 
S 4 
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3.10 Identity on the Gradient of a Product of a Vector and a Scalar 


V- (pt) = pV -t+%-Vo (3.30) 


3.11 Substantial Derivative 


Oa ane e (3.31) 


The first term on the rhs is the “local” or “unsteady” term. The second is the “convective” term. 


D 6) 6) 6) 7) 
BOBO ZO ho () +m 5e() 


“local” or “convective” 
“unsteady” (3.32) 


The rate of change D()/Dt is for two reasons: 
1. Things are changing at the point through which the element is moving (unsteady, local) 


2. The element is moving into regions with different properties. 
3.12 Differential Form of the Continuity Equation 


p(u-ri)dS + . pdV =0 (3.33) 
Ss ot Vv 


Use the Divergence Theorm {,(u- 7)dS = fff,(V - a)dV: 


SN) & ee pit) dV =0 (3.34) 


We can take the limiting case as V tends to 0 — this gives a point, where 


—+V- (pi) = (3.35) 


Figure 3.5: Taking V > 0 


Use the vector identity V - (pt) = p(V-u) +%- Vo: 


%P 4 9(V-ii) +i-Vp =0 (3.36) 
Dp oe 
ete o(V-t)= (3.37) 
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Dp Ou Ov. Ow 
= aot) = (3.38) 


Note: This means that if there is a variation of density in time and/or space, it must result in 
dilatation. One consequence is: If density is constant, then dilatation is zero. 


3.13 Momentum Conservation: Differential Form 


By arguments similar to those used with the continuity equation, the integral form of the momentum 
equation is reduced to the differential form (See textbook for derivation). Writing in terms of 
components, 


PDE made + Pfe + F’. viscous (3.39) 
Du Op 

are = dy phy F’ viscous (3.40) 
Dw Op 

PDe = az + Plz + Fz viscous (3.41) 


while f is body force per unit mass, and F' is viscous force per unit volume. In vector form, 
i dhe of 
ay ae —Vp + pf + Fryiscous (3.42) 

3.14. Energy Equation: Differential Form 


The energy equation, reduced to differential form, is: 


aa) 2 
p(e+" — 


Dt 


pP = pq—V- (pt) + of @) + Qviscous + Wviscous (3.43) 
Here, e is internal energy per unit mass, q is heat transfer rate (per unit time) per unit mass, Q 
is heat transfer per unit volume, f is body force per unit mass, W is work per unit volume. The 
equation says that the rate of change of energy per unit volume, which consists of internal energy of 
the molecules and kinetic energy of the flow, is equal to the sum of the heat transfer rate per unit 
volume, the work done by pressure on the flow (which is negative of work done by the flow against 
pressure), the work done by body forces and viscous forces, and the rate of heat transfer due to 
viscous effects, into the flow. We have not included “potential energy” due to a gravitational field 
or an electromagnetic field in the above. Gravitational potential is important when dealing with 
water flows, or with balloon flight (aerostatics). 


3.15 Simplifying the Momentum Equation 
Reynolds number: 
L datas : 
eS or Ratio of “inertial forces” to “viscous forces” (3.44) 


For typical speeds and sizes encountered in flows over aircraft, the Reynolds number is very high 
(on the order of several millions). This allows us to neglect some terms involving viscous stresses, 
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compared to other terms in the momentum equation. 


ae 
PDE 


Neglecting the viscous terms simplifies the momentum equation down to the “Euler equation.” 


SS VPs pf + F Aeeoie (3.45) 


The Euler Equation 


Di ~ 
—=-V 3.46 
PD: p+ pf (3.46) 
Where body forces can be neglected (i.e., except flows with very high radial acceleration) 
Di 
ee NT 3.47 
PD: Pp (3.47) 


3.16 Spinning Cylinder Example 


A spinning cylinder is oriented such that the rotation about its axis is 68: + 477 radians per second, 
it is in a freestream 40 km/s, with the fluid density being 1.2 kg/m°. Find the lift vector. 


L' = pU,, xT (3.48) 
7 & i 9 & 
poe =o 68 47-0 (3.49) 
0 0 40 
= 40p(47i — 68;) (3.50) 
22567 — 32647 Newtons (3.51) 


Hey, is this perpendicular to the other two vectors? Do the directions make sense? Please verify 
for yourself. 
What is the magnitude of the lift? |L’| = 2256? + 32642 =? 


3.17 Bernoulli Equation 


Derive the Bernoulli equation for steady incompressible flow, along a streamline. 
Euler equations: 


Du 


’he = —Vp (3.52) 
Ot, se Vp 
3 te VE a (3.53) 
Steady flow: 
O 
S()=0 (3.54) 
Vector identity: 
GO aie as ee (3.55) 
u = (uitvjtw rae Dy Dich 


(@-V)(@) = (ue bos ws.) @ (3.56) 
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Therefore, 


Ou Ou Ou 1 0p 
(us ee +05") meee (3.57) 
Ov | Ov Ov\ 1 Op 
(us eo, we) = ae (3.58) 
Ow Ow Ow 1 Op 
= 3.59 
(u dr Oy Oz ) poz oe) 
Multiply the first equation by dx, second by dy, and third by dz. Consider the first equation: 
Ou Ou, | Ow, lop 
up da + ae twa ie = Pr (3.60) 


Along a streamline, v/u = dy/dzx etc. So, vdx = udy, etc. 


Ou Ou Ou 1 Op 
Se aE, 3.61 
Uy ee a tur tis ee Pe ue (3.61) 
s. 1 Op 


Adding the 3 equations together, we get 


Ed aah (3.63) 
p 
Paid?) = 
dp + fa (Io ) =0 (3.64) 
p2 U3 
Pali) = 
[ rt fe fa(iaI ) =0 (3.65) 


po + FU; =pit+ Fut = const. = po (3.66) 


Chapter 4 


Potential Flow and the Laplace 
Equation 


4.1 Objective and Strategy of the Potential Flow Method 


Objective: 


Get a method to describe flow velocity fields and relate them to surface shapes consistently. 


Strategy: 


Describe the flow field as the effect of variations of one quantity. 


Vector Identity: Given a scalar function ® 
V x (V®) =0 (4.1) 


® may be ®(z, y, z, t). 
The “curl” of the gradient of a scalar function is zero. 


4.2 Velocity Potential 


We defined the “vorticity” of flow as € = V x U. This is a measure of the rotation of the flow. So 
V x U =0 means “irrotational flow.” 

Since velocity is a vector, and is a function of (2, y, z,t), this means that if the flow is irrotational, 
you can define a scalar function ®(x, y, z,t) such that 


U =Ve (4.2) 
®(x, y, z,t) is called the “velocity potential.” The definition implies: 


er 
Ox O By =e 


u (4.3) 


The condition for existence of a velocity potential is that the flow should be IRROTATIONAL. The 
speed of Mach number need not be low. 
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4.3. Incompressible Flow 


The continuity equation (conservation of mass) in differential form is 


Dp ae Dp Ou Ov Ow\ | 

Tet A(V-t)=0° or Di +p (> Dy =) = (4.4) 
Op _ Op _ Op _ Op Ou Ov Ow\ | 
ae “ae Oy GzP Oe 7 By Bz) ~ ee) 


Suppose that the changes in density anywhere in this flow are very small, and hence negligible. 


Ou Ov Ow 
0 (S++ Ge) 0 (4.6) 


If density is zero, there is no mass, and hence no flow. So the useful solution is 


0 Ov O 
ae + oF + 5: =0)| “Dilatation is zero” (4.7) 


4.4 Incompressible Potential Flow: Laplace Equation 


Incompressible: Continuity Equation is 


Ou Ov Ow 
cee i eee pe 4, 
Ox os Oy %, Oz a) 
Trrotational: 
Vx S0 (4.9) 


So the velocity can written as the gradient of a scalar Potential 


O® O® O® 
eases aes = — 4.10 
e Ox - Oy ia Oz ( ) 
Substituting, the Continuity equation becomes the Laplace Equation 
fu uv dw 
=0 Vd =0 4.11 
Ox? Oy? Oz? er ey 


So, what is the Laplace Equation, for incompressible potential flow? 
It is simply the continuity equation! Same equation describes STEADY or UNSTEADY incom- 
pressible potential flow. 

People have found several neat solutions to this equation. They have also found that it is a 
LINEAR equation, so that you can add one solution to another, and get a third solution! So we 
get the solution to flow around an airplane, as the sum of, say, 3000 little solutions to the Laplace 
equation, distributed all over the flow and airplane surface. The difficulty is in calculating the 
values of these 3000 (or whatever our computer can handle) so that their sum “obeys the boundary 
conditions” that define our flow problem. 
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4.5 Laplace Equation Approach 


For incompressible flow conditions, velocity is not large enough to cause appreciable density changes, 
so density is known — constant. Thus the unknowns are velocity and pressure. We need two 
equations. The Continuity equation is enough to solve for velocity as a function of time and space. 


e Pressure can be obtained from the Bernoulli equation, which comes from Momentum Conser- 
vation. 


e If the flow is irrotational as well, we can define a potential, and reduce the continuity equation 
to the form of the Laplace equation. 

Vv? =0 (4.12) 

O® O® O® 

—— LS =S— oS —— 

“= Oa Oy Oz 


The flow around a complete aircraft is routinely calculated from this. 


U=Ve (4.13) 


Problem: If the flow is irrotational, how can there be lift? 
Solution: We will introduce rotation, and declare the source of rotation to be “outside” the irro- 
tational flow. This strange logic works, as we will see. 


Figure 6: MGAERO-computed surface 
pressures on the Dash 8-400 
showing the aft fuselage ventral 
strakes and fin-cap fairing 


Figure 2: Panelling and computed pressures 
on the Dash 8-400, from VSAERO 


Figure 5: Grid topology for MGAERO 


computation of the Dash 8-400 www.cfdsc.ca/bulletins/14/1405.html 


Figure 4.1: Airplane Panel Code 


4.6 Building blocks 


If each of the following: ®,, ®2, 3, ®4 is a solution of the Laplace equations, then A®, + B®g + 
C®3 + D®y, is also a solution, where A, B, C, D etc. are constants. 

Thus, the solution for a complex problem can be expressed as the sum of solutions of several 
simpler problems. We will use these building blocks to analyze several flows: 


1. Uniform flow + source. 
2. Uniform flow + source + sink. 
3. As source and sink come close together, you get a doublet. 


4. Uniform flow + vortex. This looks like flow over a spinning cylinder. 
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Vortex 


Figure 4.2: Basic solutions 


We see that 1., 2. and 3. are all symmetric about the freestream direction. 4. is not. 


4.7 Boundary Conditions 


This is what specifies the details of the problem in precise mathematical terms. 


(i) Along the surface of the airfoil, the flow must be tangential to the surface: Component of 
velocity normal to the surface is zero. i.e., 


VO-n=0 U-#=0 (4.14) 


an _ 
On 


or 


0 (4.15) 


Thus, all the flow near the surface must be tangential to the surface. In other words, the 
streamline closest to the surface must be parallel to the surface. Thus another way to express 
condition (i) for 2-d problems is to say that “the surface is a streamline.” 


(ii 


Na 


Vv dy 
v _ (vy 4.16 
( U ) surface ( z) surface 


The disturbance due to an object must die away as you go far from the surface (any solution 


which says otherwise is physically unrealistic). Thus, 


As x,y,z — 00, U + Use or T — oo] +0 (4.17) 
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ii 
Figure 4.3: Far Field Boundary Condition 


Example: 
Suppose the general solution to a differential equation is ¢ = Ae’ + Be~*. What can you say 
about A or B from the above condition? 


4.8 Elementary Solutions of the Laplace Equation 


(i) Uniform flow: 


——> 
—> 


Figure 4.4: Uniform flow 


Use = . 0 (4.18) 
d=Uotf(y) ¢=const.+ g(x) (4.19) 
= fy)=comsk. oO) =i.8 (4.20) 
@ = Ucoz + const. (4.21) 


When we are using only velocity (which is the gradient of the potential), it makes no difference 
whether you have the constant: 
GO = Upet (4.22) 
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(ii) Source or sink (note: a sink is simply a negative source) 


o=const. 


Figure 4.5: Source Potential 


At the origin, there is mass flow being added or subtracted, so this is a “singular point” or 
“singularity.” There’s flow out of the plane if 2-D. Except at the origin, 


V-u=0 (4.23) 
up == ug =0 (4.24) 

r 

A 
c= (4.25) 
where . is the source strength [m?/s]. 
Og aN aN 
10¢ 
D9 7 Ue 0 @ = const. + f(r) (4.27) 
A 

= ms Inr (4.28) 


At the origin, there is mass flow being added or subtracted (“singular point”). There’s flow 
out of the plane if 2-D. 


x is a constant, related to the volume flow from/to the source/sink. 


(iii) Doublet 


Figure 4.6: Doublet Potential 
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K cosé 


o= 


On or 


(iv) Vortex 


Figure 4.7: Vortex Potential 


Let us take the “circulation” around the vortex at radius r. 


T= — pial 


T = —27rrue 
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(4.29) 


(4.30) 
(4.31) 


(4.32) 


Note that if you take the circulation at any radius, you’ll get the same value of [. The only 
vorticity is concentrated at the center: the flow is irrotational everywhere except at the center. 


OD 2, vy WOO s, 2 
ae voo Qrr 
eg 
QT 
4.9 Line Vortex 
dl 
h 
1 
P 


el 


Figure 4.8: Line Vortex 


(4.33) 


(4.34) 
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Consider a segment of a vortex filament as shown, with strength [’. Velocity induced at point P by 
the segment dl is 


ops - a (4.35) 
When @ is the angle between dl and r, 
dl x F = |dij|7| sin@ (4.36) 
Thus, velocity induced at P by an infinite, straight vortex filament is 
c= = (4.37) 


where h is the distance L” to the vortex sheet. Thus, the induced velocity drops as fh increases. 


u 


Figure 4.9: Induced Velocity 
Note: This is a general result for potential fields, taken by analogy from the result of the magnetic 
field induced by a segment dl of a conductor carrying current J, in a medium of permeability p: 


dB = ae (4.38) 
T 


4.10 Vortex Sheet 


A vortex sheet is a continuous sheet of vortices. It is used to represent a “shear layer” across which 
the velocity changes, as shown below. 


= i] 4, 
Uo a d (4.39) 
4s  {@ Vallee 


Chapter 5 


Thin Airfoil Theory 


Figure 5.1: Thin Airfoil Theory 


The lift and pitching moment produced by a thin cambered airfoil can be represented by a vortex 
sheet placed along the mean camber line of the airfoil. The issue is to find the unknown variation 
of vortex sheet strength (2). 

The boundary conditions to satisfy are: 


1) Velocity must be tangential to the vortex sheet at the vortex sheet. In other words, the vortex 
sheet is itself along a streamline. 


2) The trailing edge is a stagnation line. 


If the airfoil is represented by a vortex sheet along the camber line, then the camber line is a 
streamline of the flow. 
The velocity induced at a point x by the entire vortex sheet is 


ee aren) 
w(a) i i a (5.1) 


where € is the distance along the chord line. 

The thin airfoil assumption leads to the argument that w’(s) © w(x) where w is the component 
of velocity normal to the vortex sheet at the camber line. 
Boundary Condition: 


= 


Ux: +w'(s) =0 (5.2) 


Al 
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ws, &: = dz 
Gott = |Gre| (a - =) (5.3) 
So, note from the figure: 
> 4, we. d = dz 
eg — (Us| sin (0 — arctan =(2)) || (0 =) (5.4) 
for small angle of attack and small surface slope. 
ee ames | (2) 4 dz 
= Qn (ax pH = [ool | dx (5-9) 
dz 
we) = Us (0 - =) (5.6) 
(E)dg 
dw(x) ona — 8) (5.7) 
1 ft y(€)dé 
@=5 [TS (5.8) 
l i 
/ yg dz(x) (5.9) 
27U 5 Jo L-E dx 
cos8=1-—2zx (6 varies from 0 to 7) (5.10) 
re 
4(@) = 20 55 [46 cot 5 t oy A, sin va) (5.11) 
[cos(n — 1)” — cos(n + 1)a] = 2sin x sinnx (5.12) 
™  cosné’ F sin n@ 
[Aa = wl 
i: cos 6’ — cos 6 sin 0 ot) 
— dz 
Ao — Ay. = q— — 14 
0 » cosn0 =a Tn (5.14) 
1 [" dz 
Ap =a-—— — al 
0 a ‘ Fd (5 5) 
2 Tv 
An = — — cos nbdé (5.16) 
T JO 


If each of the following: ©, 2, ®3, ®4 is a solution of the Laplace equations, then A®, + B®:+ 
C3 + D®, is also a solution, where A, B, C, D etc. are constants. 

Thus, the solution for a complex problem can be expressed as the sum of solutions of several 
simpler problems. We will use these building blocks to analyze several flows: 


1 
2 


3 


. Uniform flow + source. 
. Uniform flow + source + sink. 
. As source and sink come close together, you get a doublet. 


. Uniform flow + vortex. This looks like flow over a spinning cylinder. 
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5.1 Results for a Symmetric Airfoil 


Symmetric airfoil: the camber line is straight. 


eS STG) 
2a Jo (x — &) 


Use the transformation: - 
f= 5 — cos 6) 


c 
c= 3 — cos 99) 


de2 5 sin od6 


Qn cos 89 — cos 00) 
7 (1+ cos 6) 
Le ree sin 0 


Consider some features of this function (5.22): 


When you go close to the leading edge, the value goes shooting up towards infinity. 


Total circulation around the airfoil is: 


r 


[reas 


= | @sinoas 
0 


Therefore, 
T = 7acU 


L! => pl Spract2: 


1 
a 5 PU oC 


cy = 27a 
dc, 


ba ey 
da : 


Pitching moment about the leading edge 


Mie =~ [Gb =—pUx fete ag 


This reduces to 


d& = Usa 


1 fn (0) sin 0 eee 
0 


a ax 9 Ta 
Mir = —doo€ ca 
The moment coefficient is ' 
= Mir =, Fe 
cea Qo 
CO 
This can also be written as ¥ 
1 
CmLE = ae 
The center of pressure is at 
/ 
jaa —-Mir Clenrr) 
ee mt + Cl 
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5.2 Cambered Airfoil 


Assume that the camber line is given by z = z(x), with z = 0 at x = 0 and x = c. The lift coefficient 
is: 


Cea It lo + - f as 69 — 1)dO (5.35) 
T Jo dx 
dc; 
ee, 
ia ot (5.36) 


Obviously, this consists of two parts: the first part is that due to angle of attack, and second due 
to camber. When the angle of attack is zero, the lift coefficient is just the second part. The angle 
of attack for zero lift is thus: 


ee pe 
Qi=0 = = ag (8 69 — 1)dO Gaya) 


We have to leave these things as integrals because the precise expression will depend on the shape 
of the camber line for the given airfoil, i.e., the function (dz/dz). 


—=C TT 
i= a = (Au eA) (5.38) 


The coefficients A; and Ag are seen from the following integrals: 


| eee ec a (5.39) 
T 0 dx 
2 Tv 
Ay.= = ae cos(n69)d8o (5.40) 
T 0 dx 
Lift and moment coefficients are: 
Woy _ 24 0, 
f= = pU2,/2 7 4 [46 cot 5 + d,An sin vs) (5.41) 
‘i A 
l= | pUdydx or cy = 2n (40 + 3 (5.42) 
0 
: 7 Ao cq Tw 
Me. = pUayxd« or aS (40 + Ay ) = (A, — Ag) (5.43) 
: 2 2 4. 4 
T 
Cm£ — — lA — Ag) (5.44) 


ep = ; c # = (At ss Ay) (5.45) 


Chapter 6 


Finite Wings 


6.1 Helmholtz Vortex Theorems 


The starting vortex is a time-rate-of-change issue. The same arguments lead to what are known as 
the Helmholtz vortex theorems for steady flow. 


1. The strength of a vortex filament is constant along its length. This leads to the next argument: 


2. A vortex filament cannot begin or end in an irrotational flow; it must extend to the boundaries 
of the fluid element, or form a closed loop. 


Figure 6.1: Helmholtz Vortex Theorems 


From the above, we can infer the vortex system generated by a finite wing moving in an otherwise 
irrotational flow. Any change in bound circulation along the wing span must be accompanied by 
the trailing of a vortex filament whose strength is equal and opposite to the said change. 


Figure 6.2: Trailing Vortex Sheet 


Thus, if the change in bound circulation is continuous, a “sheet” of vorticity, rather than discrete 
vortices, will be trailed all along the span. 


45 
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The strength of this trailing vortex sheet at any spanwise station yo is thus: 
dV 
iw) =- (5) (6.1 
1] 


The strength of each element of the trailing vortex sheet in this model, cannot change with x 
(distance downstream) unless it reaches some other flow boundary. 


y, Ty) 


Semi-Infinite Vortex Sheet 
yy) 


Visit Infinity 500 miles 
Home of the Great Starting Vortex! 


Figure 6.3: Semi Infinite Vortex Sheet 


Inboard vortex sheet 


Tip-vortices 


Figure 6.4: Vortex System of a Finite Wing 


The vortex system must extend downstream to wherever the bound circulation was first set up: 
starting vortices must have been created at the same time, and left there forever. This completes 
the vortex system of the finite wing: 


1. Bound circulation, varying along the span, and dropping to zero at each wingtip. 


2. Inboard trailing vortex sheet. 


3. Tip vortices, one from each wing tip, getting stronger downstream until all of the inboard 
vortex sheet has rolled up into the vortex on each side. 


4. Starting vortex, left behind where the bound circulation was first set up. 


In reality, the vortex sheets “rolls up” into the tip vortex as seen above. 
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6.2 Prandtl’s Lifting Line Model 
Consists of the following: 
a) A straight bound vortex I'(y), placed along the quarter-chord line of the wing. 


b) A trailing vortex distribution of strength y(y) such that 


1) == (5) oe) 


The surface boundary condition is that the normal velocity due to the wing + wake + freestream 
must be zero at the wing surface. 
In the lifting line theory, this becomes simple because: 


a) The velocity induced by the bound vortex at the bound vortex center line is zero. 


b) The wake consists of straight tip vortex + straight vortex filaments trailing back to infinity: 
semi-infinite line vortices + semi-infinite vortex sheet. 


Basis of the theory: The vortex system “induces” downward velocity everywhere inside. 


The downward velocity induced at the lifting line (wing) decreases the effective angle of attack 
seen by the wing, and hence, the lift coefficient is reduced. Also, the new “effective freestream 
velocity” seen by the wing is at an angle (the induced angle of attack) to the original freestream 
velocity. Thus the force vector is inclined backwards, and this now has a component along the drag 
direction. This is “induced drag.” 

So the problem reduces to finding the induced angle of attack, given the wing geometry, airfoil 
section characteristics, the freestream velocity, and the wing’s “geometric” angle of attack. 

In other words, 
Effective angle of attack = geometric angle of attack + induced angle of attack, where the induced 
angle of attack is negative. 

Ue = Ag + a (6.3) 


Example: Geometric angle of attack = 3 deg. Induced = -0.5 deg. Effective = 2.5 deg. 


The effective angle of attack is that corresponding to the bound circulation strength. So if we 
can find the effective angle of attack, we can find the bound circulation, and hence the lift on the 
wing. In the following, everything is written as a function of (y), the spanwise station along the 
wing. This is because the wing may be tapered (chord changes with y), twisted (angle of attack) 
or even different airfoils may be used inboard and outboard (the zero-lift angle of attack and the 
lift curve slope of the airfoil section) and anyway, the induced velocity changes with y. 


ci(y) = (ae(y) — a0(y)) (6.1) 
L'(y) = pUooT(y) (6.5) 
L! = =pU2,cly)aaly) (6.6) 
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me area hasta (6.7) 


Uso.c(y) a (y) 


The geometric angle of attack is the angle at which one sets the section including the wing twist 
if any. The induced angle of attack is the angle formed as inverse tangent of the ratio of induced 
velocity to freestream velocity. 


a; = tana! ( = ) (6.8) 


-b/2 


Figure 6.6: Spanwise Variation of Downwash 


b/2 d 
yey 
ney. = = 6.9 
w) - 4m(yo — y) 6-9) 
dV 


b/2 a dy 
Y¥ 
= —_ 4 6.10 
i An(yo — y) 
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Note on signs: 
I is positive if the wing has positive lift. ~ = dI’‘/dy is negative on the side where y is positive. 
If yo is outboard of y, yo — y is positive. The induced velocity w; at y due to the y around yg is 
downward, i.e., negative. For the point yo shown, yo — y is negative. The velocity w; at y due to 
the y around ypo is positive. 

So now our “Fundamental Equation of Lifting Line Theory” 


Me = Ag + 04 (6.11) 
can be expanded using 
1 [ UW Wi 
= x 12 
a tan (=) Us, (6.12) 
dv 
eae, /) 
wy =f (6.13) 
—»/2 47 (Yo — y) 
Lae 
yay te 
ty) + ag(y) = ag(y) + -_ / @ ud a (6.14) 
Usocly) = (y) ie a 
a 
“Fundamental Equation of Lifting Line Theory” 
dl 
bf a 2r(y) 
aol) a0) = af Get y (6.15) 
NEE Ties) Ay) 


What is unknown? Only the spanwise Bound Circulation distribution I'(y). Once we have that, we 


can find L’(y), and integrate that to get L, the lift. 


L'(y) = op ee en Gi COS Ay (6.16) 
Di(y) = Deecieely) sin ay (6.17) 

LU 1 1 
—(y) = ne — 6.18 
pi! ) tana; a; (6.18) 
Di(y) © Lay (6.19) 
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6.3 Glauert Solution to the Lifting Line Equation 


-b/2 +b/2 
y 


Figure 6.7: Glauert Solution Coordinates 


Change coordinates from y to 0. 


b 
y= 5 cos 0 (6.20) 


Assume that the unknown distribution is given by the series 


N N 
L(g) = 260. x A, sin(né) + ss Bn cnt) (6.21) 
n=1 n=1 
with the conditions 
rojeHtrniH=o (6.22) 


This works as long as the distribution does not have discontinuities — it takes only a finite number 
of terms N to represent the distribution. In practice, N = 16 is large enough. 


av d 
oe rey pg (6.23) 
89 T= BQO | Ina Jon Yo—Y 
y= > cos a= = sin 0d0 (6.24) 
at dIdo dl —2 
— — 2 

dy  d@dy (5) (saa) oe?) 

N N 
T = 20U 5, (> Ansinnd + 5 > Br cos vs) (6.26) 

n=1 n=1 


6=Oandr => I =0 (at the tips) B, =0 foralln. > T = 2bU. _, Ansinné Substituting 
into the lifting-line equation gives 


N N 


4b 
—_—_. An sin nd + ap=o(00) + nAn 
moo) Oo) ose eeet la) 2 


sin nO 


a(A9) = 


2 
sin Ao e f) 
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dc 
where mp = = at 0 = 4 corresponding to y = yo. 
a 


cosnéd6  _—_sinn®y 


Note: uses the integral = - : 
cos 0 — cos 69 sin 99 


Satisfy the relation for a(@)) at N values of #) along the span: Gives N equations for Ap, 
n=1,2,...,N. Solve. 

If we choose N points along the wing, each point is associated with one value 6. Writing this 
equation at each value of 6, we get N simultaneous equations. Solve for the N unknown coefficients 
Ay. Then 


N 
T'(8) = 2bUx. S> An sinné (6.28) 
n=1 
ue _ 28? N T 
An yo An f sin né sin 6dé 6.29 
Then 
b/2 
6.30 
oe O28 = (6.30) 
using 
- * forn=1 
| sinnésin6d0 = 4 2 (6.31) 
0 0 fornAl 
b2 
CL = Ait (5) — A,7(AR) (6.32) 


Lift coefficient depends only on the first term in the series! But you won’t get the first term right 
unless there are enough terms in the series to accurately represent the whole spanwise distribution 
of the bound circulation. 


b/2 
7 )dy 6.33 
me 0g zs | te 82) 
because 
b/2 
Di -| L'(y) sin a;(y)dy (6.34) 
—b/2 
but sina; & a; for small a;. 
on 
Cpu fia" )dy 6.35 
ve tees of 6:35) 
= 2b" A, sin n@ 9) sin 0dO 6.36 
Sg a nsinné } a;(@) sin (6.36) 
N ‘; 
= Sad (6.37) 
= sin 6 


© us = 
| sinn@sinédd =< 2 sian (6.38) 
0 0 fornAl 
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Db? pa al 9 
N A 2 
_ 2 n 
= n(AR)AZ]1+ Si () (6.40) 


Cc? ae ae 
——— — Al 
Co. = ate eee a) (6.41) 
Ch Ch 
= 1 = 6.42 
AO) ARE Oe) 
6.4 Elliptical Lift Distribution 
i Ge Ayn(AR) 
— — A 
D6, ; N A,\? (6.43) 
A2n(AR) [1+ a ($) 
L 
Max — = A, =0forn> 2. 
dD; 
Dy \2 
= [= 2bU.A1 sin 8 = 2bU oo AVI — cos? 8 = 2bU or 4/1 — (=) 
2y\? : : 
T=[94/1- oe Equation to an ellipse! (6.44) 
For this I’ distribution, downwash is 
W (60) = 5 ; constant along the span! (6.45) 
—W Cr 
i= = A 
a v.. > a(AR) (6.46) 
Ch 
= =1 6.47 


Another simplification is possible if the distribution is symmetric: no yaw or roll, and you’re 
not carrying a spare engine on one wing, for instance. 


T(6) = T(r -8@) (6.48) 
An = 0 for even values of n. So, 
N 
DT = 2bU x0 S~ Aon—1sin(2n — 1)8 (6.49) 
nt 


Thus with just 4 points, you can get a series with 7 terms — up to A7. 


Chapter 7 


Viscous Flow 


7.1 Viscous Shear Stress 


Ou Ov 

wes “lt oe) 
Ov Ow 
“(St a) 


OD si0e) 
a sae FS 


Qq 

xR 

8 
| 


(7.1) 
(7.2) 


(7.3) 


The direction is indicated by the second subscript. The units of stress are Force/Area, such as 


N/m?, psi, or psf. 


These can be written in compact form. Here u is used to represent any of the velocity com- 
ponents, and zx is used to represent any of the spatial coordinates. i and j representing x,y,z in 


turn. 


—— Ou; + Ou; 
na B Ox; Ox; 
Normal Stresses due to Viscosity 


6) ) re) O 
Ora = aioe + ( aes aa ) 


Ox Oy Oz 
Ci SH 2 os +A Ou + Oe + au 
me Poy Ox Oy Oz 


Stokes’ Hypothesis: 


So the normal stress due to viscosity add up to zero. 
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(7.4) 


CHAPTER 7. VISCOUS FLOW 54 


7.2 2D Navier-Stokes Equations For Incompressible Flow 


Conservation of mass: 3 3 
ub wv 
—+—=0 7.9 
Ox - Oy oo 


Conservation of z-momentum: 


Ou Ou 1LdOp | Ou Ou 
Conservation of y-momentum: 
Ov Ov 10p au Ou 
| Sy a 
“Oa "By p Oy (3 xt) aay) 


7.3 Exact Solution to the Steady 2D Navier-Stokes Equations For 
Incompressible Flow 


Plane Parallel Steady Flow Between Two Infinite Plates: Couette Flow 


y, V h 


Figure 7.1: Couette Flow 


From the problem, there cannot be any change in u along x: So du/dz is zero. Also, v = 0. So 
u is u(y). 


1 dp du 
a) ee T12 
pdx Yay? ey) 
1dpy? 
u(y) = a +Ay+B (7.13) 
pdx 2 


2 
= Aye (7.14) 
a 
Boundary conditions: u=v=0aty=0,uw=U andv=0O0aty=h. 


U 


y, V 


Figure 7.2: Flow Between Driven Plates 
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y dp Uy 
= | wl 
u(y) Dy 7 ts (7.15) 
If dp/dx = 0, 
iy) = Linear velocity profile (7.16) 
Shear stress s F ” 
u v 
ry = = 7.17 
Skin Friction Coefficient ‘ 
ee mM 
= = 7.18 
spU2 phU ‘Re ee) 
puh 
Re = Reynolds number based on U and h = ae (7.19) 
Case 2: Flow between two stationary plates, driven by pressure gradient 
yy, V 
x 
Figure 7.3: Flow Between Stationary Plates 
y dp 
SS ph 2 
u(y) See (y —h) (7.20) 


Note: Sign of dp/dz must be opposite to that of u. 


Hagen-Poiseuille Flow: Steady axisymmetric, fully-developed incompressible lam- 
inar flow through a constant-diameter pipe 


du du 10 FG) seu | 

"5, +9, t oae =” bi a m8 ea) 
dv dv 10 12(r3") a 
ug tea, 4 cae oe oe eg) 


u-momentum equation reduces to 


du 
10p | wf (“S) 


poz rr adr 


(7.24) 
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Integrating, 
d 1 
u dp Ea C 
dr 2p dz r 


For flow to remain finite at axis (r = 0), C =0. Also, u = 0 at the wall, r= R. 


1 dp 


a (R?—r?) Parabolic Velocity Profile 
pdz 


u(r) = 
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(7.25) 


(7.26) 


Chapter 8 


Boundary Layer 


8.1 Boundary Layer Equations 


Boundary Layer Approach 


Newton’s 2nd Law 


| 


Conservation Equations in integral form 


| 


Convert to differential form 


iE 
Simplify by use of OM analysis 
Boundary Layer Approximations 


At very high Reynolds number (~ 100, 000+) based on distance along the surface in the freestream 
direction, 


1. The spatial rate at which properties change across a boundary layer is very large, compared 
to the rate at which things change along the flow direction. In other words, derivatives with 
respect to y are much higher than derivatives with respect to 2. 


d( )/dy >> d( )/dax (8.1) 
2. Order of magnitude analysis shows that the v-momentum equation reduces to 


Op 
EES 8.2 
ied (8.2) 


In other words, static pressure is constant across a boundary layer. 


u-momentum equation 


2 2 
Ou Ou ae aE | (8.3) 


v | = 
Ox Oy pdx p |x? 5 Oy? 
For high Reynolds number flows, our order of magnitude analysis says that the u-momentum equa- 


tion may be approximated as 
Ou Ou | 10p wu 


1 T CN. A 
Ox "Oy pox poy? ee 
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v-momentum equation 


becomes ap se oe 
dy ; 
au a ~ =0 (8.7) 
3a yt abe” aot 88) 
) 
a ~0 (8.9) 


a) The variation of pressure across the boundary layer is small. Thus, we can assume that pressure 
varies only with respect to x. The pressure field may be computed at the edge of the boundary 
layer from inviscid flow theories, without any knowledge of the boundary layer characteristics. 


b) The v-component of velocity within the boundary layer is of the order of (U.6/c) where 6 is the 
boundary layer thickness. 


c) The boundary layer thickness 6 is very small, and for laminar flows it varies inversely as the 
square root of Reynolds number pUc/u. 


d) We will assume the airfoil surface to be flat, and use a Cartesian coordinate system, neglecting 
surface curvature. This is because the boundary layer thickness 6 is very small compared to the 
airfoil surface radius of curvature. 


8.2 Boundary Layer Features 


Velocity profile 


Figure 8.1: Boundary Layer Features 
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1. Boundary Layer Thickness 6: 
Defined as the y—location where u/ue reaches 0.99%, that is the u-velocity becomes 99% of 
the edge velocity. 


2. Displacement Thickness 6*: 
This is a measure of the outward displacement of the streamlines from the solid surface as a 
result of the reduced u-velocity within the boundary layer. This quantity is defined as 


ee i c = | dy (8.10) 
0 Pete 
3. Momentum Thickness 6: 


This is a measure of the momentum loss within the boundary layer as a result of the reduced 
velocities within the boundary layer. 


o= fe )1- 4] dy (8.11) 


Pete 


4. Shape Factor H: 
This quantity is defined as the ratio 6*/0. 
For laminar flows H varies between 2 and 3. It is 3.7 near separation point. Thus excessively 
large values of H (above 3) indicate that the boundary layer is about to separate. In turbulent 
flows, H varies between 1.5 and 2. 


8.3. Wall Shear Stress and Drag 


Surface Shear Stress: 
The shear stress at the wall can be found from definition of shear stress 


Ou 
ie fee eta deta 8.12 
wall HM ( ay ) aii ( ) 


Skin Friction Coefficient cy: 
Tw 


cr = cies 
pPertle 
Skin Friction Drag D: 


Shear stress may be numerically integrated over the entire solid surface to give the skin friction 
drag force along the x-axis 


(8.13) 


| pele (8.14) 
Entire Surface 


Skin Friction Drag Coefficient Cg: 
The drag force is usually non-dimensionalized by the freestream dynamic pressure times the chord 
of the airfoil c, giving the skin friction drag coefficient along the x-axis, Cg. 


(8.15) 
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8.4 Thwaites’ Integral method for Laminar Incompressible Bound- 
ary Layers 


This is an empirical method based on the observation that most laminar boundary layers obey the 
following relationship (Ref: Thwaites, B., Incompressible Aerodynamics, Clarendon Press, Oxford, 
1960): 

Ue (07) 6? due 


=A-B Ll 
vy dx vy dx (8-16) 


Thwaites recommends A = 0.45 and B = 6 as the best empirical fit. 
The above equation may be analytically integrated yielding 


Oday $f? ae 
62 = a uedx + e252 (8.17) 


For blunt bodies such as airfoils, the edge velocity ue is zero at x = 0, the stagnation point. For 
sharp nosed geometries such as a flat plate, the momentum thickness q is zero at the leading edge. 
In these cases, the term in the square bracket vanishes. 

The integral may be evaluated, at least numerically when we is known. 

After @ is found, the following relations are used to compute the shape factor H and the shear 
stress at the wall 7,. For 0 <A < 0.1 


H= 9.61 =3.75\-+4 5.24)? (8.18) 

For —0.1 <A <0 aie 
H = 2.472 + —____ 8.19 
+ 0107+ ery 

where is 
Ue 

yeas 8.20 
vy dx ( ) 
hie “0 + 0.09)°:82 (8.21) 


Despite the empiricism in the above formulas, Thwaites’ integral method is considered to be the best 
of integral boundary layer methods. 


8.5 Transition 


Consider the 2-D unsteady, incompressible viscous flow past an airfoil at a sufficiently high Reynolds 
number 


U 
Re = —* & 300, 000 (8.22) 
m 


We find that the boundary layer over both the upper and lower surfaces may be divided into three 
regions. They are (1) laminar region, (2) transitional region, and (3) turbulent region. 
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Transition strip: rough strip taped to airfoil 
close to max thickness location, to trigger distur- 
bances that cause quick transition to turbulence. 


Figure 8.2: Transition 


8.6 Turbulence Features 


The turbulent region has unsteady flow, although the fluctuations are small and occur about some 
“steady” mean flow levels. The velocity profile is fuller than in the laminar flow case. Skin friction 
is higher than in laminar flow. 


u, Vv, p 


Mean Value 


Time 


_ Turbulent 


1 : i ; 
| Potential Core: | Transition Zone | 


h ages Laminar Geen ' Fully Developed 


x - Turbulence 


Figure 8.3: Turbulence Features 


CHAPTER 8. BOUNDARY LAYER 62 


8.7 Incompressible Boundary Layer on Flat Plate at zero angle of 
attack: Laminar vs. Turbulent 


Laminar Turbulent 

5.02 0.372 

= Ros d= pee 
_ 0.6642 _ 0.664 
Re®5 °f Re®5 

1.328 0.074 

Cr= Cr = = 

f Re®5 f Re®-? 


8.8 Incompressible Turbulent Skin Friction Correlation At Large 
Reynolds Numbers 


Schultz-Grunow equation, modified by Kulfan, Boeing Co. 
Cr. = 0.295 (low Hes) 7s. 10° & Re, < 10" (8.23) 


0.0050 


& 
a 
n 


*log(Re x * Cfi) + 1.7 


Cfi 


0.0040 


IC fi = 0.295 * (Log(Re x)) > 


0.0030 


0.0020 


Incompressible Local Skin Friction ,Cfi 


0.0010 


0.0000, 95 407 


© Theory: Modified Schultz / Grunow 
= Theory: Karmen - Schoenerr 


Reynolds Number, Rex 10° 10° 


Courtesy, Dr. R. Kulfan, Boeing Co. 


Figure 8.4: Comparison of Incompressible Local Skin Friction Predictions 


Chapter 9 


Some Practical Results 


9.1 Airfoil Drag vs. Angle of Attack 


Note: When viscous drag and/or flow separation are present, airfoils do have drag that changes 
with angle of attack. In other words, C'pg of an aircraft is not completely independent of angle of 
attack, over a large range of angle of attack. This has nothing to do with tip vortices etc. So when 
we must use an averaged value for C'p9 when we use lifting line theory and say that for a whole 
aircraft (note nomenclature). 


QO 


Figure 9.1: Airfoil Drag Coefficient at Angle of Attack 


CG 
Ll 
Cp = Co t+ iRye (9.1) 


The more usual and general way, is to represent the drag as 


Cp =Cpo+ KC} (9.2) 


K will now be MARE 
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Lmax 


Figure 9.2: Aircraft Drag Polar 


9.2 Relation between Boundary Layer Profile and Skin Friction 
Drag 


A positive pressure gradient dp/dz would cause the boundary layer velocity profile to become 
“unhealthy” — develop an inflection point. 
Beyond that, 3 things can happen: 


1. Laminar separation 
2. Transition to turbulence 


3. Turbulent separation 


_~ Turbulent 


aoe Laminar 
QNNN NNN NN? * 


YNNNN 


Figure 9.3: Positive pressure gradient would make boundary layer velocity profile “unhealthy” 
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Ou 
TET i 9.3 
wall ‘ (x ) ai) 03) 


Inflection 
point 


Stable: Unstable Turbulent 
No inflection point 


Figure 9.4: Transition to Turbulence 


Turbulent boundary layer has high velocities close to surface, large du/dy — high skin friction drag 
(bad). However, very stable profile delays separation (good!). 


9.3. Separation 


3-D Separation Line 


Nal _ 


Stagnation dy =0 
wall 


Figure 9.5: Flow Separation 
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9.4 Separation Control 


Figure 9.6: Separation Control: surface blowing (upper) and cavities (lower) 


9.5 Moving Wall 


Wall Moving 


Fixed Wall Downstream 


Wall Moving Downstream 


et ag 


Figure 9.7: Moving Wall Effect 
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9.6 Boundary Layer Suction 


Figure 9.8: Boundary Layer Suction 


9.7 Vortex Interactions & Ground Effect 


Technique: Think of the velocity induced at the axis of each vortex by every other vortex. A solid 
wall is treated as an “image plane” with an imaginary mirror image of the vortex on the other side. 


1 
errs, 
- 
° 1 wy 
wgers SS 

ole * 

1 f ores 
Cr 2 ee 

3 ao o - 


ee 
: 
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Figure 9.9: Vortex Interactions and Ground Effect 
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9.8 Potential Flow Panel Code with Interacting Boundary Layer 


Calculation 


“XFOIL” by Prof. M. Drela  http://web.mit.edu/drela/Public/web/xfoil/ 


E 184 Re= 100000 


Althaus Experiment. 


Ais 0.020 cy 9.025 


oy. SN pe 
Se : 
—-— 


Figure 9.10: Panel Code With Boundary Layer 


9.9 Drag Coefficients 


Shape Effects on Drag 


The shape of an object has a very great 
effect on the amount of drag. 


Cd=1.28 Cd=1.14 


Flat Plate 


a 
a. Cd= .295 
Bullet 


Flow | 


Cd=.07 to 5 
Sphere 
Cd = ene Cd= .045 
p V“A/2 Airfoil 


A= frontal area All objects have the same frontal area. 


Figure 9.11: Shape Effects on Drag 
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Table 9.1: Drag Coefficients for bluff bodies 


Rough sphere (Re = 10°) 0.40 
Smooth sphere (Re = 10°) 0.10 
Hollow semi-sphere opposite stream 1.42 
Hollow semi-sphere facing stream 0.38 
Hollow semi-cylinder opposite stream 1.20 
Hollow semi-cylinder facing stream 2.30 
Squared flat plate at 90° 1.17 
Long flat plate at 90° 1.98 


Open wheel, rotating, h/D = 0.28 0.58 


Table 9.2: Drag Coefficients for streamlined bodies 


Laminar flat plate (Re = 10°) 0.001 
Turbulent flat plate (Re = 10°) 0.005 


Airfoil section, minimum 0.006 
Airfoil section, at stall 0.025 
2-element airfoil 0.025 
4-element airfoil 0.05 
Subsonic aircraft wing, minimum 0.05 
Subsonic aircraft wing, at stall 0.16 
Subsonic aircraft wing, minimum 0.005 
Subsonic aircraft wing, at stall 0.09 
Aircraft wing (supersonic) n.a. 


Table 9.3: Vehicle Drag Coefficients 


Subsonic transport aircraft 0.012 
Supersonic fighter, M=2.5 0.016 


Airship 0.020-0.025 
Helicopter download 0.4-1.2 
Sports car 0.3-0.4 
Economy car 0.4-0.5 
Pickup truck 0.5 
Tractor-trailer, with fairings 0.6-0.7 
Tractor-trailer 0.7-0.9 
Trailer alone 0.9 

Racing car 0.65-1.10 


The most famous source on drag data is: 
“Fluid-Dynamic Drag: Practical Information on Aerodynamic Drag and Hydrodynamic Resistance” 
by Sighard F. Hoerner, Dr.-Ing.-habil. 


